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The experimental signatures of nonstandard neutrino interactions are shown to be equivalent to 
CPT violation in neutrino oscillations. This result leads to a correspondence in the study of these 
two descriptions that can be used to constrain the relevant parameters of one formalism by using 
the available bounds on parameters of the other. The correspondence is illustrated and explicitly 
used to determine first bounds on previously unexplored parameters in both formalisms. 


I. INTRODUCTION 


Neutrino oscillations have been experimentally con¬ 
firmed using accelerator, atmospheric, reactor, and solar 
neutrinos [1]. These observations indicate the remark¬ 
able interpretation of massive neutrinos. The minimal 
model extending the Standard Model (SM) to accommo¬ 
date neutrino masses has been verified with great preci¬ 
sion in all the experiments above [1]. Even though to date 
the absolute mass scale of neutrinos remains unknown, 
the interferometric nature of neutrino oscillations has al¬ 
lowed us to measure the effects of their minute masses 
in the form mass-squared differences. This unexpected 
behavior within the context of the SM has motivated 
the search for unconventional effects that could manifest 
their minuscule effects in neutrino oscillations. 

These so-called exotic scenarios beyond the SM have 
become very active fields and their study gave rise to an 
interesting incidental program of experiments originally 
designed to measure the parameters of the three-neutrino 
massive model. Some of these ideas include the search for 
nonstandard neutrino interactions (NSI) [2, 3], Lorentz 
and CPT violation in neutrinos [4], long-range interac¬ 
tions [5], large extra dimensions [6], and sterile neutrinos 
[7]. Regarding the first two exotic scenarios mentioned 
above, the physical motivations are completely indepen¬ 
dent; nonetheless, the observable effects that they would 
produce in experiments can be related. This means that 
although they remain as independent modihcations of the 
conventional physics, the experimental constraints of one 
can help to constrain the other. 

In this work, a correspondence between matter NSI 
and CPT violation in neutrinos is presented, which allows 
relating the parameters in these two formalisms. This pa¬ 
per is organized as follows. Matter NSI and their effects 
in neutrino oscillations are described in Sec. II, while the 
corresponding effects of CPT violation are presented in 
Sec. III. In Sec. IV the correspondence between mat¬ 
ter NSI and CPT violation is discussed and the explicit 
application of a perturbative method is illustrated. The 
proposed correspondence is used in sections V and VI to 
relate and determine new bounds on the parameters of 
both formalisms. Sec. VII describes a brief discussion on 
the sensitivity prospects for future experiments. 


II. NONSTANDARD INTERACTIONS 

Neutrino interactions at low energies can be effectively 
described by four-fermion vertices containing two neu¬ 
trino states. In particular, nonstandard interactions of 
neutrinos with quarks and electrons in a material medium 
can be described by the effective Lagrangian [2, 3] 

Unsi = -H h.c. (1) 

where Gp is the Fermi constant, P = L,R, and the 
strength between neutrino states v of flavors a and (3 and 
the L-handed (i?-handed) components of the fermions / 

f f f f 

and /' is parametrized by the factors ). For 

f ^ f, the NSI produce charged-current type effects that 
modify the neutrino production and detection. These ef¬ 
fects are controlled by NSI parameters denoted and 
referring to the source and detector, respectively. Re¬ 
markably, these parameters can trigger neutrino flavor 
change even in the limit of no propagation through the 
so-called zero-distance effect [8] and several proposals for 
the experimental study of these parameters have been 
considered [9]. In the present work, we are interested on 
NSI parameters that modify the neutrino propagation 
through matter, which occurs in the case / = /' in the 
form of NSI neutral-current type effects. 

In vacuum, the effective hamiltonian describing neu¬ 
trino oscillations can be written in terms of the mass- 
squared differences Am^i, Arri^i, and the neutrino en¬ 
ergy E in the form 

1 /O 0 0 \ 

iJo = —C/ OAmil 0 [/t, (2) 

yo 0 Am^iJ 

where U is the PMNS matrix that parametrizes the 
mixing between flavor and mass eigenstates in terms 
of four constant parameters: three mixing angles and 
one CP phase [10]; Majorana phases could also be in¬ 
cluded, although they are unobservable in neutrino os¬ 
cillations. For neutrino experiments on Earth, in par¬ 
ticular long-baseline beam experiments, the propagation 
through matter can play a crucial role [2]. Forward scat¬ 
tering of neutrinos with the particles that they encounter 
as they propagate in a dense medium can be modelled as 
a constant contribution to the hamiltonian in the form 
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H = Hq-\- Hm, where 


I 0 0 \ 

Hm = V^Gprie 0 0 0 

\o 0 oy 


( 3 ) 


This matter potential arises from charged-current inter¬ 
actions with electrons in the medium of number density 
rie- Neutral-currents modify the three neutrino flavors by 
the same amount; being proportional to the identity in 
flavor space, these interactions produce no effects in the 
oscillations of three neutrinos so they can be disregarded. 

In a similar fashion, the vector component of the NSI 
(1) for f = f modifies neutrino propagation through 
matter, in which case the combination 
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controls the observable effects. In the last expression, 
the number density of the relevant fermion of type / is 
denoted by n/ and we have considered that neutrinos 
propagate through an unpolarized medium. The physical 
effect of the parameters is a modification of the neu¬ 
trino dispersion relation similar to the matter effect (3) 
affecting all flavor components independently. The cor¬ 
responding effective hamiltonian takes the general form 

( pTn _m pm \ 

^ee ^e/j. '^er \ 

C . (5) 

_m* _m* ~m I 

^er ^ (IT ^rr/ 

where the flavor structure guarantees the hermiticity of 
the hamiltonian. For the propagation of antineutrinos, 
the corresponding NSI hamiltonian is obtained by the re¬ 
placement —)■ —-\/2Gi?ne£™^. Notice that 

due to the absence of antimatter in the medium, spurious 
CP- and CPT-violating effects can appear when neutri¬ 
nos propagate in matter [11]. 

The full hamiltonian describing neutrino oscillations 
through a dense medium in the presence of NSI can be 
written as 


H = Ho + Hm + H^si- (6) 

The study of the parameters in the last term is ad¬ 
dressed by the diagonalization of the full hamiltonian (6). 
Given the mixed energy dependence of the different terms 
in this hamiltonian, the mixing angles will be energy- 
dependent. This phenomenon is well known from the 
study of solar neutrinos, in which case the matter poten¬ 
tial is given in terms of the electron density in the solar 
core. Moreover, the eigenvalues of the hamiltonian will 
exhibit an unconventional dependence of the neutrino en¬ 
ergy that will differ from the vacuum behavior that makes 
the oscillation phase proportional to E~^. Different ap¬ 
proaches have been implemented for the search of matter 
NSI using atmospheric [12], beam [13], reactor [14], and 
solar [15] experiments. Distinct methodologies have been 


explored for the direct study of the NSI parameters in 
neutrino oscillations, including approximation techniques 
that can be implemented in some experimental configura¬ 
tions. In particular, given the remarkable success of the 
model of three massive neutrinos, NSI are usually consid¬ 
ered as sub-leading effects that could modify the conven¬ 
tional description of neutrino oscillations [9, 12-15]. To 
date, most of the constraints on the NSI parameters £™^ 
have been obtained from phenomenological studies [16]. 
Two direct experimental studies have been performed by 
the MINOS [17] and Super-Kamiokande [18] collabora¬ 
tions. Hereafter, we only consider bounds obtained from 
experimental studies. 


III. CPT VIOLATION 

Neutrino oscillations are natural interferometers that 
offer great sensitivity to search for new physics, including 
deviations from exact Lorentz invariance. The potential 
breakdown of one of the most fundamental symmetries 
of modern physics has been mostly motivated by string- 
theory scenarios [19]. Systematic searches of deviations 
from Lorentz symmetry have been implemented using 
the so-called Standard-Model Extension (SME) [20, 21]. 
These modern tests of Lorentz invariance span a wide 
range of fields, whose experimental results are tabulated 
in Ref. [22]. In the neutrino sector [4], different tech¬ 
niques have been developed for weak decays [23-25] and 
astrophysical neutrinos [26]; nevertheless, neutrino oscil¬ 
lations using short- and long-baseline experiments have 
shown to be very sensitive probes of unique forms of 
Lorentz violation that lead to neutrino mixing [27, 28]. 

There exists a subset of Lorentz-violating operators 
that also break CPT invariance in the fermion sector. 
Restricting attention only to renormalizable Dirac cou¬ 
plings in the theory, these CPT-odd operators are written 
in the form [4, 29] 

/ICPT- = - + ^.C., (7) 

with the flavor indices taking the values a, (3 = e, /i, r 
and where the coefficients and control the im¬ 
pact of vector and pseudo-vector couplings, respectively. 
It is important to emphasize that coordinate invariance 
is preserved and only particle Lorentz invariance is bro¬ 
ken [30]. We remark in passing that Majorana couplings 
can also be included in the lagrangian (7), which can 
mix neutrinos and antineutrinos [4, 28]. These Majorana 
couplings produce direction-dependent effects and have 
been experimentally constrained using beam [31], reac¬ 
tor [32], and double-beta-decay experiments [24]. From 
the lagrangian (7), the observable effects on left-handed 
neutrinos are controlled by the combinations 

{aL)al3 = (a + (8) 

which are constant hermitian matrices in flavor space 
that modify the conventional vacuum hamiltonian 
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through an observer-independent structure of the form 
(®i)a/ 3 -PA [4, 29]. The corresponding hamiltonian for 
antineutrinos is obtained after the coefficients (8) are 
substituted by the right-handed coefficients = 

(a — 6)^^ = —The spacetime index A = 0,1, 2, 3 
explicitly exhibits the potential breakdown of rotational 
invariance, where = (l;p) depends on the neutrino 
direction of propagation p. Nonetheless, in the present 
work we are only interested in the isotropic component 
iaL)%. 

The explicit form of the CPT-violating modification of 
the neutrino hamiltonian can be written as 

/(aL)L {aLfer 

HcpT- = («L)°r 

\{aL)Tr iaL)% {aLfrr 

where the components of completely character¬ 

ize independent deviations from CPT invariance. The 
full hamiltonian describing oscillations of CPT-violating 
neutrinos can be written as 

H = Ho + Hm + Hcpt-, ( 10 ) 

where we have included the matter hamiltonian for com¬ 
pleteness. The diagonalization of this hamiltonian can 
then be used to study the coefficients ( 0 ^)°^ in oscillation 
experiments. For this purpose, perturbative methods 
have been developed for short- and long-baseline experi¬ 
ments [27, 28]. Many experimental studies have been per¬ 
formed to constraint these coefficients by Double Chooz 
[33], IceCube [34], LSND [35], MiniBooNE [36, 37], MI¬ 
NOS [38-40], and Super-Kamiokande [41]. Similarly, sen¬ 
sitivity studies have been performed for the future JUNO 
experiment [42]. The hamiltonian (10) has also been used 
to implement alternative models for neutrino oscillations 
[43]. 

IV. NSI-CPT-VIOLATION CORRESPONDENCE 

Comparing the effects at the hamiltonian level of NSI 
(5) and CPT violation (9) we can directly make the pa¬ 
rameter correspondence 

icLL)^P ( 11 ) 

This expression establishes an equivalence between the 
observable effects introduced by the matter NSI param¬ 
eters and the coefficients for CPT violation (ol)®^. 
The modified neutrino oscillation probabilities in both 
scenarios will take the same form when describing CPT- 
violating neutrinos modelled by the theory (7) or when 
introducing matter NSI described by the theory (I); 
the relation (11) shows the connection between the two 
formalisms. This result implies that the experimental 
bounds obtained in one formalism could be translated 
into bounds for the other. 

It is important to emphasize that even though the re¬ 
lation (II) allows us to relate the parameters of matter 


NSI and CPT violation in neutrino oscillations, the un¬ 
derlying physics controlled by the corresponding param¬ 
eters remains inequivalent. In fact, matter NSI require 
that neutrino propagate through matter to be observable, 
whereas the effects of CPT violation are independent of 
the medium. In both cases, there exists an effective re¬ 
fractive index that modifies the neutrino dispersion re¬ 
lation. For matter NSI, this index of refraction arises 
due to the unconventional interactions with electrons and 
quarks in the medium, which act like a condensate that 
alters the neutrino propagation. On the other hand, the 
index of refraction in the CPT-violating scenario corre¬ 
sponds to the existence of an intrinsic background field 
that isotropically permeates the vacuum. The nature of 
this an other background fields has been extensively stud¬ 
ied for theories with Lorentz invariance violation. For an 
updated review, see Ref. [44] and references therein. 

Now that we have established a direct relationship be¬ 
tween the matter NSI parameters and the coefficients 
that control CPT violation, we can borrow the meth¬ 
ods from one framework to apply it to the other. Cau¬ 
tion is necessary for a correct use of the relation (II) 
because there might exist situations in which the rela¬ 
tionship between parameters is invalid due to the nature 
of the analysis to extract limits on the parameters. For 
instance, some bounds on NSI have been obtained assum¬ 
ing a particular matter composition so the corresponding 
electron density must be treated with care. Conversely, 
some bounds on coefficients for CPT violation have been 
obtained using neutrinos that propagate in different me¬ 
dia; hence, again the correct interpretation of the matter 
density must be taken into account. 

Many of the constraints on CPT-violating neutrinos 
have used perturbative methods, treating CPT viola¬ 
tion as a sub-leading effect over the conventional mass- 
driven oscillations, an approach also implemented in the 
study of matter NSI. For instance, for long-baseline ex¬ 
periments the far detector is located at a distance from 
the source that maximizes the conventional oscillation ef¬ 
fects. A power-series expansion of the oscillation proba¬ 
bility in the form + 

... allows identifying the different terms in the series as a 
function of the sub-leading physical process added to the 
conventional description. As expected, the first term in 
the series is the conventional neutrino oscillation proba¬ 
bility driven by neutrino mass-squared differences. Using 
the results in Ref. [28], this conventional probability can 
be written as where the zeroth-order 

oscillation amplitude for neutrinos that propagate a dis¬ 
tance L is given by 

si°i=J2Uo.kU*0ke-^^'‘^. (12) 

k 

In this expression, the mixing matrix Uak and the eigen¬ 
values Ek of the unperturbed problem are necessary. If 
conventional matter effects are negligible, then the mix¬ 
ing matrix and the eigenvalues are simply given in terms 
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of the of the mixing angles and mass-squared differences 
of the vacuum hamiltonian (2). 

The second term in the series corresponds to the first 
modification due to the unconventional physics, which 
according to (11) could be either NSI or CPT violation. 
For isotropic CPT violation, this term has the explicit 
form [28] 


p(i) 

up- 


Ua. 


= 2LIm 





7(5 


(13) 


where the sum over flavor indices shows that different 
components of the hamiltonian (9) can be studied with 
this method. The complex weighting factors 
are defined in terms of the unperturbed mixing parame¬ 
ters [28]. The probability (13) is linear in the coefficient 
controlling CPT violation and arises due to the inter¬ 
ference between the conventional hamiltonian and the 
perturbative correction (9). Making use of the relation 
(11), the above expression for the probability can then be 
used to study the components of the NSI hamiltonian (5). 
The overall factor L in the probability (13) shows that 
the minute effects of the parameters ( 0 ^)°^ or can 
be enhanced by a long propagation distance, as expected 
from an interferometric measurement. 

The second-order term in the perturbative expan¬ 
sion of the oscillation probability corresponds to two 
quadratic contributions [28]. The first of these effects 
is analog to the linear probability (13) arising from the 
interference between CPT violation and the mass-driven 
hamiltonian. The second quadratic effect is only due to 
CPT violation, leading to oscillations even if the con¬ 
ventional effects are negligible. This is the leading-order 
term in the perturbative series for experimental config¬ 
urations for which the zeroth-order amplitude is negligi¬ 
ble. This condition applies to experiments having short 
baseline compared to the mass-driven oscillation length, 
either because the propagation distance is too short or 
because the neutrinos used have very high energy [27]. 

In addition to perturbative methods, the exact treat¬ 
ment of the components of the coefficient ( 0 ^)°^ can also 
be implemented. In fact, for experiments studying neu¬ 
trinos over a large range of energies and baselines, such 
as atmospheric neutrinos, perturbative methods can be¬ 
come unapplicable. Exact diagonalization methods have 
been available in the literature since the realization of 
the importance of matter effects. Many of these meth¬ 
ods have been designed for long-baseline experiments, in 
which case neutrinos travel for several hundreds of kilo¬ 
meters through rock [45]. For completeness, a very com¬ 
pact method for the exact diagonalization of an arbitrary 
3x3 hermitian hamiltonian as well as the general form 
of the mixing matrix are presented in Appendix A. 


V. BOUNDS ON NSI FROM CPT VIOLATION 

Most of the searches for CPT violation in neutrinos 
have focused on coefficients that generate sidereal vari¬ 


ations of the oscillation probability. Only a few ex¬ 
periments have studied spectral distortions generated 
by isotropic coefficients for CPT violation. The Mini- 
BooNE collaboration obtained the 2a bound |(aL)g^| < 
4.2 X 10“^° GeV [36, 37] and the Double Chooz collabo¬ 
ration determined the limit \{aLY^r\ < x I0“^° GeV 
at the 95% G.L. [33, 46]. In both cases, neutrinos prop¬ 
agate only a few hundred meters through the Earth’s 
crust, where the matter density is about 1 gr/cm^ [47]. 
Given the relatively short baseline of these two exper¬ 
iments, their sensitivity is limited and the listed limits 
correspond to the bounds on NSI parameters 


|e™|<I.lxI0^ I <2.0x10% (14) 


These bounds are orders of magnitude weaker than the 
phenomenological limits already known in the literature 
[16]. Nevertheless, it should be emphasized that the val¬ 
ues (14) are the first bounds on the magnitude of the 
and e™ components directly obtained from an experi¬ 
mental analysis. 

In a recent study, the Super-Kamiokande collaboration 
determined the most stringent upper bounds on real an 
imaginary parts of the off-diagonal components of the 
isotropic coefficient for CPT violation at the 95% C.L. 
[41] 


Re(aL)°^ < 1.8 X 10-23 GeV, 
lia{aL)% < 1.8 X 10-23 GeV, 

Re(aL)°^ < 4.1 x 10-23 GeV, 

Im(aL)°^ < 2.8 X 10-23 GeV, 

Re(aL)°^ < 6.5 X 10-24 GeV, 
lm{aL)lr < 5.1 X 10-24 GeV. (15) 


It is tempting to simply apply the relation (II) to trans¬ 
late the bounds on ( 0 ^)°^ into however, given the 
multiple directions and creation points of atmospheric 
neutrinos used for the analysis, the matter potential ex¬ 
perienced by these neutrinos is far from unique. A mod¬ 
est estimate can be obtained by using an appropriate 
average density for all the neutrino events. The data 
analysis found that up-going events are the most sensi¬ 
tive sample to the unconventional effects studied, which 
is expected due to the large travel path of these neutri¬ 
nos. Since 90% of up-going neutrinos propagate through 
matter with an average density of 3.4 gr/cm3 or more 
[47], we can take this value as a moderate approximation 
for the average density of matter used to determine the 
bounds (15). We remark that this is a conservative es¬ 
timate because a significant fraction of the events used 
in the analysis traverse denser layers of the Earth that 
would lead to even more stringent bounds on the NSI 
parameters. The relation (II) then gives the following 
upper bounds on the real and imaginary parts of the NSI 
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FIG. 1: Individual existing and new experimental bounds on 
the matter NSI parameters The two rows for the cases 
a ^ P correspond to real and imaginary parts, respectively. 
The estimated sensitivity of the future DUNE experiment is 


also shown. 



parameters 



Re£™ < 1.4 X 

1 

o 

Ime™ < 1.4 X 10-3 

Ree™ < 3.2 x 

I 

o 

Im£™ < 2.2 X 10-3 

Re£™ < 5.1 x 

fXT 

10-2, 

Ime"i < 4.0 X 10-2 

f-lT 


It should be noticed that the values (16) include the first 
bounds on the imaginary parts of these off-diagonal pa¬ 
rameters. Until now, only the magnitude of these param¬ 
eters has been considered. Additionally, these bounds 
are the first on the and e™ components directly ob¬ 
tained from an experimental analysis. Even considering 
indirect phenomenological bounds on these components, 
the values (16) exhibit a slight refinement on |e™ |, and 
improvements by more than a factor five on \e^^\ and al¬ 
most one order of magnitude on |e™ |. Fig. 1 illustrates 
existing bounds as well as the improvements obtained in 
this work by considering that only one of the parameters 
at a time is nonzero. 


VI. BOUNDS ON CPT VIOLATION FROM NSI 

Systematic experimental searches of matter NSI have 
been only performed to study the /j,t sector. Using ac¬ 
celerator neutrinos, the MINOS collaboration [17] im¬ 
plemented a two-flavor parametrization of NSI and using 
both neutrino and antineutrino data determined a bound 
on the real part of at the 90% C.L. given by 

-0.20 < Re e™ < 0.07. (17) 

Reaching a maximum depth of about 10 km, more than 
81% of the 735-km path of the neutrino beam of the MI¬ 
NOS experiment goes through rock of density 2.6 gr/cm^ 
[47]. Thus, the two-sided bound (17) translates into a 
limit on the coefficient for CPT violation 

-1.95 X 10-23 GeV < Re(aL)° < 6.83 x lO-^^ GeV. 

(18) 

The remaining 19% of the beam propagates through less 
dense matter; however, we have only kept the most con¬ 
servative value. Notice that the upper bound is slightly 
weaker than the limit shown in (15); nevertheless, the 
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FIG. 2: Individual existing and new experimental bounds on 
the coefficients for GPT violation (ar)))^ in GeV. The two 
rows for the cases a ^ P correspond to real and imaginary 
parts, respectively. The estimated sensitivity of the future 
DUNE experiment is also shown. 


result from Super-Kamiokande is one sided, whereas the 
result (18) includes the first lower bound on Re(ai)°^. 

Similarly, using atmospheric neutrinos and a two- 
flavor approximation for the characterization of NSI, the 
Super-Kamiokande collaboration determined the follow¬ 
ing bounds in the fj,T sector at the 90% C.L. [18] 

|Re£™ I < 1.1 X 10-2, 

I < 4.9x10-2. (19) 

Following the approach presented in Sec. V, atmospheric 
neutrinos propagate through matter of different density; 
therefore, an appropriate application of the correspon¬ 
dence in Eq. (11) requires the use of a reasonable value 
for an average density. In the same conservative approx¬ 
imation used in Sec. V, the two-sided bounds (19) trans¬ 
late into the following limits on the corresponding coef¬ 
ficients for CPT violation 

|Re(ai)°J < 1.4 X 10-24 GeV, 

- (.aL)°r\ < 6.3 x 10-2^ GeV. (20) 

In addition of being two-sided bounds, the above result 
provides a slight improvement on Re(ai)3^ and the first 
limit involving diagonal components. Following the stan¬ 
dard method for estimating attainable sensitivities, we 
consider one coefficient at a time to write the following 
bounds on the relevant diagonal components |(aL)°^| < 
6.3 X 10-24 GeV and \{aL)°r \ < 6-3 x 10-24 GeV. Fig. 2 
illustrates existing bounds as well as the improvements 
obtained in this work by considering that only one of the 
parameters at a time is nonzero. 

VII. FUTURE PROSPECTS 

Since the neutrino propagation distance enhances the 
minute effects of e™ and (ul)® , several current and fu¬ 
ture experiments, such as the Deep Underground Neu¬ 
trino Experiment (DUNE) [48], NuMI Off-Axis Ve Ap¬ 
pearance experiment (NOi^A) [49], and the Precision Ice- 
Gube Next Generation Upgrade (PINGU) [50] offer ex¬ 
citing prospects for the study of matter NSI and GPT 
violation. 
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As an illustration, we can determine an estimate on 
the sensitivity of DUNE. A comprehensive sensitivity as¬ 
sessment can be done using GLoBES [51]; however, here 
we simply consider a crude estimate in the form of a 10% 
deviation from the conventional oscillation probability. 
DUNE is a multipurpose experiment designed for a rich 
physics program including the study of long-baseline neu¬ 
trino oscillations, supernova neutrinos, and atmospheric 
neutrinos. As a long-baseline experiment, improvements 
on the current bounds are expected due to the 1300 km. 
that separate the source and the far detector. Sensitiv¬ 
ities to the matter-NSI parameters are presented in 
Fig. 1, while the corresponding sensitivities to the co¬ 
efficients for CPT violation ( 0 ^)°^ are presented in Fig. 
2. In both cases, the potential bounds are illustrated by 
considering that only one of the parameters at a time is 
nonzero. 

For matter NSI, we find that the bounds from CPT 
violation obtained in this work filled most of the previ¬ 
ously unexplored parameter space with great sensitivity. 
Nonetheless, DUNE could improve the bounds on the 
electron sector. For CPT violation, bounds on all the 
coefficients in the electron sector could be improved by 
DUNE as well as the imaginary part of ( 01 ,)°,.. 


VIII. SUMMARY 

Neutrino oscillations offer a remarkable sensitivity to 
search for new physics. Two particular descriptions of 
unconventional physics have been explored and shown 
to produce similar modifications. Despite the com¬ 
pletely different theoretical motivations and underlying 
physics, nonstandard matter interactions and CPT vi¬ 
olation, parametrized by the matrices and ( 0 ^)°^, 
respectively, lead to the same observable effects. Since 
the data analysis only makes use of the oscillation prob¬ 
abilities, from an experimental point of view the results 
from one framework are equivalent to the other. This 
equivalence is explicitly presented in Eq. (11), which has 
been used to relate the current bounds the parameters 
on both formalisms. 

Under mild assumptions, new limits have been deter¬ 
mined in these sets of parameters derived from a system¬ 
atic experimental analysis. Despite the stringent phe¬ 
nomenological bounds on |, \s'^\, and in Eq. 

(16) these limits have been improved and bounds on 
the imaginary parts of these parameters have been con¬ 
strained for the first time. Similarly, an increased sensi¬ 
tivity and a lower bound on the coefficient (ul)®,- is found 
in Eq. (20) as well as the first bound on the combination 
of diagonal components [( 0 ^)°^ — {aL)%\- 

These results demonstrate the functionality of the cor¬ 
respondence between formalisms presented in this work. 
Future bounds on coefficients for CPT violation in neu¬ 
trino oscillations can be used to report bounds on matter 
NSI parameters, and vice versa. In the case of a positive 
signal, a distinction between matter NSI and CPT vio¬ 


lation can be made because the former can only occur 
for neutrinos that propagate through a dense medium, 
whereas the observable effects of the latter can occur in 
vacuum. 

The vast experimental program in neutrino oscillations 
offers exciting opportunities for the study of neutrinos 
as well as the search of unconventional physics. In the 
present work we have established a direct correspondence 
between two communities in the hunt for new physics 
using neutrinos. 
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Appendix A: Exact diagonalization 

We begin with the 3x3 hermitian hamiltonian in flavor 
basis, whose entries will be denoted by Hap. Denoting 
the three eigenenergies by Ej and the identity matrix by 
1 , the eigenvalue equation det{H — Ejl) =0 leads to the 
cubic equation that can be nicely written in terms of the 
hamiltonian invariants 


E^ MtEj MttEj -\- Md = 0 , (-^ 1 ) 


where 


Mt = -TtH, 

Mtt = i(TrEr)2-iTr(J?2), 

Md = —deiH. (A2) 


The matrix invariants are functions of the entries of the 
hamiltonian, which can be used to implement Cardano’s 
method for solving the cubic equation (Al). By defining 
the angles 


arccos 


'2Ml - QMtMtt + 27Md \ 


2(M2 - dMTTf/^ ) ^ 


(A3) 

with ti = 0, ^2 = —tz = 27r, the three solutions of the 
cubic equation (Al) can be written as 


E,= 


2(M| - dMTTf/'^ cos6»j -b Mt 


(A4) 


Once the eigenenergies have been found, we can use the 
eigenvalue equation to construct the eigenstates, which 
will constitute the columns of the mixing matrix. The 
nine entries of the mixing matrix can be explicitly written 
as 


Uej — 


N, 




Nj \aj 


U, 


TJ 


l«jl 
N, 1 / 3,1 ’ 

(A5) 
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where we have defined the functions 

I3j — H ^/j,r-^rej 

7i = H^r{Hee- Ej) - H^eHer, (A6) 

and the normalization factors 

N, = (|a,/3,f + |a,7,f + l/3,7,f (A7) 

Now that we have explicit forms for the eigenenergies 


(A4) and the elements of the mixing matrix (A5), we 
can use the general form of the oscillation probability 

= d„;3-4;^Re(C/AC/^,[/„fcC/^fc) sin2(lA,fct) 

j>k 

+2^Im(C/AC/^,C/„fcC/^fc) sin(A,fet) (A8) 

j>k 

to write the probabilities of interest, where the eigenvalue 
differences are A^/j = Ej — Ek- 
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